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Abstract 

Band diagrams often pass through special k-vector points with “irrational” fractional coordinates such as 
some Brillouin zone vertices and centers of edges. Use of “standard” primitive cells as defined by Setyawan 
and Curtarolo [Comp. Mater. Sci. 49, 299] is convenient for band diagram and effective mass calcu¬ 
lations because the definitions and coordinates of special k-vector points including “irrational” fractional 
points are provided, but their basis vectors are taken differently and/or special k-vector point definitions 
differ from the crystallographic convention in many cases. On the other hand, the Bilbao Crystallographic 
Server [Bulg. Chem. Commun. 43, 183] defines special k-vector points based on crystallographic con¬ 
vention; however, symbols for “irrational” special k-vector points are not dehned. Obviously there is a 
need for a crystallographic convention-friendly scheme to describe such “irrational” special k-vector points. 
Therefore, we derived a procedure to uniquely choose basis vectors for all space groups from the crys¬ 
tal symmetry, basis vector lengths, and inter axial angles, and implemented this algorithm in the spglib 
code I http://spglib.sourceforge.net/ . This paper outlines an algorithm that is relatively straightforward 
to implement in software to bridge between crystallographic conventional cells and “standard” primitive 
cells by retaking basis vectors. The procedure is expected to ease systematic electronic and phonon band 
calculations. 


1. Introduction 

An algorithm for efficient and systematic genera¬ 
tion of electronic and phonon band diagrams would 
be useful not only in detailed investigation of cer¬ 
tain compounds but also in combinatorial or “high- 
throughput” calculations UM- Band diagrams and 
effective carrier masses are typically obtained by 
sampling k-points along a path spanning the first 
Brillouin zone (BZ). Although the shape of a BZ 
is independent on how the basis vectors are cho¬ 
sen, the fractional coordinates of vertices and cen¬ 
ters of edges and faces, hereafter denoted as spe¬ 
cial k-vector points, are not. Therefore, rules to 
determine basis vectors that give the same special 
k-vector point coordinates for a given form of a BZ 
would be useful. The choice of basis vectors, def¬ 
initions of special k-vector points, and suggested 
band diagram paths are provided in detail for “stan¬ 
dard” primitive cells as defined by Setyawan and 
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Curtarolo [3]. However, the “standard” cell differs 
from the crystallographic conventional cell in many 
cases, and for this reason we use the term “stan¬ 
dard” with quotation marks in this article when 
referring to Setyawan and Curtarolo’s definition. 
Consequently, there is demand for a simple proce¬ 
dure, preferably using transformation matrices for 
retaking basis vectors, to convert a crystallographic 
conventional cell into a “standard” cell. We note 
that the AFLOW code Q can automatically gen¬ 
erate “standard” cells; however, details of the algo¬ 
rithm to derive a “standard” cell is not disclosed. 

The objective of this article is to review how the 
conventional and “standard” definitions are differ¬ 
ent and to outline an algorithm to convert a crys¬ 
tallographic conventional unit cell of the standard 
setting into a “standard” primitive cell. This algo¬ 
rithm is comprehensive and relatively easy to imple¬ 
ment in software. Section [5] provides a summary of 
definitions and notations used in this paper. Sec¬ 
tion [3] introduces the crystallographic convention 
on basis vector choice. Section |4] outlines the dif- 
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ferences between crystallographic and “standard” 
conventional cells for each Bravais lattice and gives 
procedures and transformation matrices to convert 
the former to the latter. Section [5] outlines how 
“standard” primitive basis vectors are expressed in 
Cartesian coordinates. The definitions are based 
on Setyawan and Curtarolo but re-written in a 
form more convenient for implementation in soft¬ 
ware. Section M shows criteria that determine the 
topology of the BZ. Section [7] brings up the neces¬ 
sity of handling special k-vector points with “irra¬ 
tional” coordinates and discusses issues arising in 
comparing special k-vector point symbols between 
those in the Bilbao CrystalloCTaphic Server 0 , 0 , @] 
and Setyawan and Curtarolo (3| . Section[8]is a sum¬ 
mary of this paper. 

2. Symbols and terminologies 

Definitions of symbols are summarized in Table 
[T] Symbols for basis vectors are written in this 
article as ai, a. 2 , and as in line with the Fifth edi¬ 
tion of the International Tables of Crystallography 
A (ITA) 01 and the Third edition of the Interna¬ 
tional Tables of Crystallography B (ITB) 0. The 
“standard” conventional and primitive cells are de¬ 
fined according to Setyawan and Curtarolo 0. In 
essence, reciprocal space basis vectors of “standard” 
primitive cells penetrate the center of Bragg planes 
belonging to the first BZ and direct space basis vec¬ 
tor lengths are chosen as a' <b' < c' when there is 
ambiguity. Primes are added for “standard” cells 
and the subscript P indicates a primitive cell. All 
direct space basis vectors are column vectors and 
reciprocal space basis vectors are row vectors. The 
direct and reciprocal space basis vectors are related 
by • aj = Sij , • a^- = Sij , ap^ • apj = 6ij 

{i,j G {1,2,3}) using the Kronecker delta symbol 
Sij. Lengths of basis vectors and interaxial angles 
are collectively referred to as lattice constants. 

The symbols of the Bravais lattices and the sec¬ 
tions where topics are discussed in this article with 
respect to each Bravais lattice are summarized in 
Table H 

A change of basis matrix M relates a choice of 
basis vectors to another choice of basis vectors. M 
is a 3 X 3 matrix and is used as 

(a'i,a2>a3) = (ai>a2,a3)M. (1) 

A transformation matrix, P, transforms a conven¬ 
tional unit cell to one of the primitive cells by 

(api, ap 2 , apg) = (a}, a 2 , a 3 )P. (2) 


With the exception of rhombohedral cells, the 
change of basis M conserves the number of atoms 
in the unit cell, ria,, because det(Ai’) = 1. On the 
other hand, det(P) = 1 for cells without centring, 
det(P) = 1/2 for base- and body-centred cells and 
det(P) = 1/4 for face-centred cells. The number of 
atoms in the primitive cell is ria det(P). 

The unit cell is a parallelepiped determined by 
the basis vectors 0 . A point in the unit cell x is 
represented by a column vector of fractional values 
with respect to basis vector lengths: 

= {xi,X 2 ,X 3 ) where 0 < Xi < 1. (3) 

The position of an atom at a point x in the unit 
cell is given by 

X = ^Xiai = (ai,a2,a3)a;. (4) 

i 

X is invariant with change of basis. However, as 

(ai,a2,a3)a; = {a[,a2,a'2)M~^x, (5) 

the vector representation of the atom varies with 
change of basis M according to 

x' = M~^x -\-1, (6) 

where t is the lattice translation vector necessary to 
fulfill the condition m so that the point is in the 
new unit cell. 

Irrational special k-vector points are defined as 
special k-vector points that always have irrational 
fractional coordinates irrespective of basis vector 
lengths except in virtual crystals. Some vertices 
of the BZ are located at irrational special k-vector 
points in body-centred tetragonal, face-, body-, and 
base-centred orthorhombic, rhombohedral, mono¬ 
clinic, base-centred monoclinic, and triclinic Bra¬ 
vais lattices. 

3. Convention on basis vector choice 

There is a crystallographic convention on which 
specific choice of basis vectors should be used out 
of an inhnite number of choices, and this is out¬ 
lined in the ITA 0 and ITB 0. Here, the Bra¬ 
vais lattice is determined, symmetry operations are 
identified, and then the origin is decided. A total 
of 530 choices, distin^ished by Hermann-Mauguin 
and Hall symbols 0 , 0 , of 230 space group types 
are listed in the ITA and ITB. Multiple choices are 
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Table 1: List of symbols used in this study. 


Definition 

Basis vectors 

Basis vector lengths 

Interaxial angles 

Crystallographic conventional cell 

{ai, a2, as} 

{a,b,c} 

{a,/3,7} 

Reciprocal crystallographic conventional cell 


kc\ 


“standard” conventional cell 

{a'i,a^,a(j} 

{a',6',c'} 

W,/3',7'} 

Reciprocal “standard” conventional cell 

{a'r,a-,a-} 

{K,K,K} 


“standard” primitive cell 

{api, ap2, apg} 

{Op,&p,Cpj 

K,/3',7'} 

Reciprocal “standard” primitive cell 

{apg, ap2, apg} 

{kpa, fcp^, fcp^} 

{fcp„, kpg, kp^} 


Table 2: Symbols of the Bravais lattices and sections of this article where they are discussed. 


Crystal system 

Centring 

Symbol 

Sections 

Triclinic 

Primitive 

aP 

[ssiiiTiini] 

Monoclinic 

Primitive 

mP 

17^0110 


Base-centred 

mS 


Orthorhombic 

Primitive 

oP 

[TT1IT31I6T] 


Base-centred 

oS 



Body-centred 

ol 

ITTIOIIOTIIO 


Face-centred 

oF 

1770115.2.4116.31 

Tetragonal 

Primitive 

tP 

177177110] 


Body-centred 

tl 

[74117717771721 

Hexagonal 

Primitive 

hP 

[77177177 

Rhombohedral 

* 

hR 

177177177 

Cubic 

Primitive 

cP 

177177177 


Body-centred 

cl 



Face-centred 

cF 

1771771777177 


^Primitive or triple hexagonal. 


provided for many space group types. For instance, 
space group type Pnn2 (No. 34) can be Pn2n 
or P2nn depending on how the basis vectors are 
chosen, and space group type P4,/ncc (No. 130) 
has two choices with different origins. In general, 
changing the basis may change the centring in base- 
centred cells but does not in body- or face-centred 
cells; one exception is that body-centred monoclinic 
cells can be converted into base-centred cells by us¬ 
ing a different choice. The standard choice for mon¬ 
oclinic cells always has unique axis b and [3 > 90° 
and base-centred cells are always C-centred. Base- 
centred orthorhombic cells are C-centred with the 
exception of four space group numbers 38 to 41; 
the conventional setting of these four space group 
types is A-centred such that the point group is mm2 
and the proper two-fold rotation axis is along the 


c-axis. The triple hexagonal unit cell for rhombo- 
hedral Bravais lattices could be obverse or reverse; 
the current version of the International Tables takes 
the obverse setting as the standard setting although 
the reverse setting was the standard until 1952 [7|. 

The ITA by itself does not uniquely define the 
crystallographic conventional cell in a number of 
situations. For instance, the order of basis vector 
lengths is not specihed, which becomes problem¬ 
atic in many monoclinic and orthorhombic space 
group types. One standardization algorithm to 
uniquely determine the conventional cell is pro¬ 
posed by Parthe et al. [13 They require the cell to 
be right-handed and have the following constraints: 
the Niggli reduced cell is used for triclinic struc¬ 
tures, “best” cell with unique axis b and /? > 90° 
[HI is used for monoclinic structures, that is, the 


3 


























cell where a and c are as short as possible while con¬ 
forming to the standard space group symbol, a < b 
or a < b < c for orthorhombic structures, and ob¬ 
verse triple hexagonal cell for rhombohedral struc¬ 
tures. Additionally, origin 2 is used for centrosym- 
metric structures. The space group type with the 
smallest index for the relevant screw axis is taken. 
For instance, one can choose to define space group 
number 80 as but, as this space group type 
can also be defined as /4i, the former is not ac¬ 
cepted as a choice of the space group type and the 
latter is adopted because the screw axis index is 
smaller. Parthe et al. 0 provides further rules 
to reduce the values of internal coordinates that 
we will not discuss in this paper. The above pro¬ 
cedure unambiguously defines the basis vectors of 
the crystallographic conventional cell but the in¬ 
ternal coordinates are not necessarily determined. 
One example is CsCl with space grou p ty pe Pmim 
as mentioned in Parthe and Gelato [l^. Wyckoff 
positions la and 16, which both have the same site 
symmetry mSm, must be occupied by Cs or Cl, but 
the above requirement does not determine whether 
Cs should occupy la or 16 sites. This ambiguity 
still remains after the standardization procedure of 
Parthe et al. [l^ However, the band diagram is 
not affected by the choice of origin since this is a 
physical property in reciprocal space that is invari¬ 
ant under translation in direct space. The conven¬ 
tional cells employed in this work conform to the 
conventional coordinate system in Table 2.1.2.1 of 
the ITA [^. The basis vectors cannot be uniquely 
defined from Table 2.1.2.1 in many space group 
types, and therefore the following additional rules 
are imposed. We take the first choice that appears 
in Table Al.4.2.7 of the ITB Q as the standard 
choice. This means that origin 1 is always used and 
the space group type with the smallest index for 
the relevant screw axis is automatically taken. All 
other choices in a particular space group type can 
be converted into the first choice by change of basis 
and/or translation of the origin. Table 9.3.7.1 of the 
ITA provides definitions of conventional cells based 
on lattice constants; however, the resulting space 
group choice is not necessarily the first one in Ta¬ 
ble Al.4.2.7 of the ITB Q, therefore the definition 
in Table 9.3.7.1 of the ITA is not adopted. Parthe 
et al. 0 instead takes the symmetry center at the 
origin, that is, adopt origin choice 2, for centrosym- 
metric structures. Origin choice 2 is never the first 
choice that appears in Table Al.4.2.7 of the ITB Q. 
Parthe et al. [10| also requires that “orthorhombic 


structures where the space group symbols do not 
prescribe a particular labelling of any of the three 
unit-cell axes have to be described with a unit cell 
where a < 6 < c”; however, this condition is inap¬ 
propriate in some orthorhombic space group types 
as discussed below. Conventional cells are defined 
as above in the spglib code [l3| versions 1.7.3 and 
higher. 


3.1. Orthorhombic lattices 

Some space group types have ambiguity on how 
to take basis vectors. An obvious example is space 
group type P222 (No. 16) where the six permuta¬ 
tions of a, 6, and c do not change the space group 
type. 

The conditions to resolve ambiguity in basis vec¬ 
tor lengths of orthorhombic space groups, with or 
without centring, other than Pbca (No. 61) and 
Ibca (No. 73) are a < b, and if necessary, addition¬ 
ally 6 < c. For example, only the former is nec¬ 
essary in Pcc2 (No. 27) and both is necessary in 
Fmmm (No. 69). The imposed condition for Pbca 
and Ibca is to make a the shortest. Two choices 
for Pbca exist: Pbca (adopted choice) and Pcab. If 
Pbca is a < 6 < c then Pcab is a < c < 6, and if 
Pcab is a < 6 < c then Pbca is a < c < 6. There¬ 
fore, a < b < c cannot be enforced if Pbca is to 
be used as the adopted choice. Similar logic ap¬ 
plies to Ibca (No. 73). The point groups that are 
allowed for an orthorhombic lattice are 222, mm2, 
and mmm. As the proper two-fold rotation axis 
is taken along the c-axis in the point group mm2, 
Amm2, Abm2, Ama2, and A6a2 (numbers 38 to 
41) become A-centred. 


3.2. Monoclinic lattices 

The definition according to Parthe and Celato 


111 is adopted and is summarized as follows. The 
6-axis is taken as the unique axis and /3 > 90°, and 
the basis vectors corresponding to a and c are taken 
as short as possible. Base centring, if it exists, is 
always C-centred. The condition a < c is imposed 
if the c-glide symmetry element identihable by the 
“c” in the Hermann-Mauguin symbol is absent. 


3.3. Triclinic lattices 

The crystallographic conventional cell for a tri¬ 
clinic system is the Niggli reduced cell that can be 
robustly obtained using the algorithm by Kfivy and 
Cruber Q and its modified version that addresses 
numerical errors arising in finite precision floating 
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point algebra by Grosse-Kunstleve et al. [l^ Parthe 
et al. [l0| also employs the Niggli-reduced cell. 


3 . 4 . Cubic, tetragonal, hexagonal, and rhombohe- 
dral lattices 

For cubic, tetragonal, and hexagonal lattices, the 
usual conditions are adopted when choosing basis 
vectors in this article, that are, 

a = b = c and a = /3 = 7 = 90° for cubic, 
a = b and a = (3 = "f = 90° for tetragonal, 
a = b,a = (3 = 90° and 7 = 120° for hexagonal, 

respectively. For primitive rhombohedral cells, the 
condition is 


easily from the space group number. The cell is 
triclinic if the space group number is 1 or 2 , mono¬ 
clinic if 3 to 15, orthorhombic if 16 to 74, tetragonal 
if 75 to 142, hexagonal or rhombohedral if 143 to 
194, and cubic otherwise. The first character of the 
Hermann-Mauguin symbol gives the centring (P, 
A, C, I, F, or R) and distinguishes hexagonal and 
rhombohedral cells for space group numbers 143 to 
194. 

4 . 1 . cP, cF, cl, tP, tl, and hP 

The crystallographic and “standard” Q conven¬ 
tional cells are the same for these Bravais lattices; 
therefore, the change of basis matrix M is the iden¬ 
tity matrix. 


a = b = c and a = (3 = ^. 

However, the conventional rhombohedral cell is de¬ 
fined on a hexagonal lattice. The obverse setting 
of the conventional rhombohedral cell can be re¬ 
covered from the primitive rhombohedral cell using 
the transformation matrix of Eq. © shown 

in Sec. 14.21 

4. Deriving basis vectors of the “standard” 
cell 

The crystallographic conventional unit cell as de¬ 
fined in Secs. [XT] to [Til is chosen as the starting 
point except in triclinic cells as discussed in section 
14.71 Computer programs typically represent basis 
vectors in Cartesian coordinates and store internal 
coordinates of atoms using fractional coordinates or 
Cartesian coordinates. In this case, converting be¬ 
tween fractional and Cartesian coordinates is easily 
done by a linear transformation such as Eq. @. 
Standardization of basis vectors is carried out by 
choosing different basis vectors, and the positions of 
atoms are invariant in Cartesian coordinates during 
retaking of basis vectors. The volume of the unit 
cell is kept constant with the exception of rhombo¬ 
hedral cells. 

We will not discuss cases when the basis vectors 
cannot be uniquely defined. This can happen in 
virtually designed crystals; one example is an or¬ 
thorhombic cell with space group type 1222 and 
a = b = c. The basis vectors cannot be uniquely 
defined unless a, b, and c are all different from each 
other in this case. 

The first step of the outlined algorithm is identifi¬ 
cation of the Bravais lattice, which can be deduced 


4.2. hR 

The rhombohedral crystallographic conventional 
cell is obverse triple hexagonal while the “standard” 
conventional cell is primitive rhombohedral[^. The 
change of basis matrix and the inverse are 



333 


/l 0 l\ 

M = 

112 

333 

,M-i = 

I 1 1 


U - -} 

\3 3 3/ 


^0 i 1 ) 


respectively. The determinant of M is 1/3, which 
means that the “standard” conventional cell vol¬ 
ume is one-third of the crystallographic conven¬ 
tional cell. One simple procedure to remove dupli¬ 
cate atoms in the “standard” cell is to apply Eq. 
and remove duplicate atoms. 

4-3. oP, oF, and ol 

The order of basis vector lengths is not always 
a < 6 < c in a non-base centred crystallographic 
conventional cell whereas the order in a “standard” 
conventional cell is always a' < 6 ' < c' Q. The 
change of basis matrix M and the inverse are given 
in Table [3] 

/./. oS 

A base-centred orthorhombic crystallographic 
conventional cell is generally C-centred but some 
space group types with point group type mm 2 , 
namely space group types Amm2, Abm2, Ama2, 
and Aba2 (38 to 41) are A-centred. The order of ba¬ 
sis vector lengths is not unique in some space group 
types. In contrast, a “standard” conventional cell 
is always C-centred and a' < b' regardless of sym¬ 
metry [Sj] . The change of basis matrix M and the 
inverse are given in Table |4l 
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Table 3: Change of basis matrix M for non-base 
centred orthorhombic cells. 


Table 4: Change of basis matrix M for base centred 
orthorhombic cells. 


Condition 


a < b < c 


a < c < b 


b < a < c 


b < c < a 


c < a < b 


c < b < a 


M 

JTT 

0 1 
^0 0 
A 0 
0 0 
Vo 1 
/o 1 
1 0 
Vo 0 
/o 0 
1 0 
Vo 1 
/o 1 
0 0 
Vl 0 

/o 0 

0 I 

yi 0 



M -1 

/l 0 o\ 
0 10 
Vo 0 ij 
/l 0 o\ 
0 0 1 
Vo 1 oj 
(o 1 o\ 
10 0 
Vo 0 Ij 
(o 1 o\ 
0 0 1 
Vl 0 o) 
/o 0 l\ 
10 0 
Vo 1 oJ 
(o 0 l\ 
0 I 0 

yi 0 0/ 


4-5. mP 

As outlined in Parthe and Gelato 0, a crystal¬ 
lographic conventional cell has unique axis b and 
/3 > 90°. The additional constraint of a < c is 
always imposed on space group types without c- 
glide symmetry. The order of a and c is not al¬ 
ways the same when there is c-glide symmetry be¬ 
cause this symmetry element distinguishes ai and 
as. On the other hand, a “standard” conventional 
cell has unique axis a, lattice constants b' < c', and 
a' < 90° [3|. The change of basis matrix M and 
the inverse are given in Table [SJ 


4-6. mS 

A crystallographic conventional cell is C-centred 
with unique axis b and /3 > 90°, which is required 
in a “best” cell as defined by Parthe and Gelato 


llj. In contrast, a “standard” conventional cell 


is C-centred with unique axis a and a' < 90° 

We believe the condition h' < c' in the original pa¬ 
per by Setyawan and Curtarolo cannot be enforced 
because the basis vector length of in the plane 
with centring is not necessarily shorter than the ba¬ 
sis vector length of ag without centring. The three 


Condition 


C-centred, a <b 


C-centred, b < a 


A-centred, b < c 


A-centred, c<h 


M 

/l 0 o\ 
0 10 
Vo 0 l) 
/o 1 o\ 
10 0 
Vo 0 1 ) 
/o 0 l\ 
10 0 
Vo 1 0/ 
/o 0 l\ 
0 10 
yi 0 0/ 


M -1 

/l 0 o\ 
0 10 
yo 0 i) 

(o 1 o\ 
10 0 
yo 0 i) 

/o 1 o\ 
0 0 1 
Vl 0 0 / 

(o 0 l\ 
0 10 
yi 0 0/ 


Table 5: Change of basis matrix M for monoclinic 
cells. 


Condition 


a < c 


c < a 


M 

/o i o\ 
10 0 
yo 0 ij 
/o 0 l\ 
loo 

yo 1 0/ 


/o 1 o\ 
loo 
yo 0 ij 
/o 1 o\ 
0 0 1 
yi 0 0/ 


basis vectors are distinguished by symmetry; one of 
the basis vectors comprising the plane with centring 
forms a right angle with the unique axis while the 
other does not. The change of basis matrix and its 
inverse are 



/o 

I 

0\ 

(0 

1 

of 

M = 

1 

0 

0 ,M-^ = 

I 

0 



lo 

0 

1/ 

lo 

0 

1/ 


The reason why only one change of basis matrix is 
sufficient is provided in [Appendix A 


4.7. aP 

The crystallographic conventional triclinic cell is 
the Niggli reduced cell [13, IS Ell that is deter¬ 
mined with direct space lattice constants. The or¬ 
der of basis vector lengths are a < b < c, and a, /3, 
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and 7 are all smaller than (Type I) or larger than 
(Type II) 90°. On the other hand, the “standard” 
conventional cell is defined based on reciprocal basis 
vectors. The requirement that reciprocal basis vec¬ 
tors must penetrate centers of faces of the first BZ 
means that the reciprocal “standard” conventional 
cell must be a Niggli reduced cell in reciprocal space 
with axes relabeled. In other words, the set of basis 
vector lengths of the reciprocal “standard” conven¬ 
tional cell must match that of the Niggli reduced 
cell in reciprocal space. Moreover, the interaxial an¬ 
gles of the reciprocal “standard” conventional cell, 
k'^, k'p, and fc(,, must be either all larger than or 
smaller than 90° and k'^ is always the one closest to 
90° 0. 

The conversion from crystallographic to “stan¬ 
dard” cells is done in the following three steps. 
First, the Niggli reduced cell in reciprocal space is 
obtained and its reciprocal basis vectors, that is, 
basis vectors in direct space, are used as ai, a 2 , 
and aa instead of the basis vectors of the crystal¬ 
lographic conventional cell. Next, basis vectors are 
relabeled such that k^ is closest to 90° by using a 
change of basis matrix Mi that transforms to an 
intermediate form: (a",a 2 ,a 3 ) = (ai, a 2 , a 3 ) 7 Vfi . 
Finally, another change of basis matrix M2 is used 
as (a'a, a2, a'^) = (a", a^', a!^)M2 to invert basis vec¬ 
tors as necessary such that ka, k^, and k^ become 
all larger than or smaller than 90°. In other words, 

= (ai>a2,a3)MiM2. (9) 

The change of basis matrix for the second step, Afi, 
is given in Table |6l The change of basis matrix for 
the third step, M2, depends on the interaxial angles 
of the reciprocal unit cell of the intermediate form 
fc", kp, and fc" and is given in Table [71 


Table 6 : Change of basis matrix Mi for triclinic 
cells. 


Condition 


Ml 


Mi 


-1 


k' closest to 90° 


k'a closest to 90° 


k!y closest to 90° 



Table 7: Change of basis matrix M 2 for triclinic 
cells. 


Condition 


/c" < 90°, k'^ < 90°, k” < 90° or 
/c" > 90°, > 90°, /c" > 90° 

/c" < 90°, fc" > 90°, k” > 90° or 
/c" > 90°,fc(j' < 90°,/e" < 90° 

/c" > 90°, fc" < 90°, /c" > 90° or 
/c" < 90°, fc" > 90°,/c" < 90° 

/c" > 90°, fc" > 90°, fc" < 90° or 
fc" < 90°,fc;j' < 90°, fc" > 90° 


M2 = Mi'i 

/l 0 o\ 
0 10 
Vo 0 ij 
/l 0 o\ 
0 I 0 
Vo 0 Ij 
(1 0 o\ 
0 10 
Vo 0 Ij 
(1 0 o\ 
0 I 0 

yo 0 V 


5. ’’Standard” basis vectors in Cartesian co¬ 
ordinates 


5.1. ’’Standard” conventional cell 


“Standard” basis vectors can be given in Carte¬ 
sian coordinates once the “standard” lattice con¬ 
stants are obtained. Representation of points in a 
unit cell as given in Eq. (jS)) is invariant against rigid 
rotation of the corresponding crystal. The follow¬ 
ing are how “standard” basis vectors, both conven¬ 
tional and primitive, are represented in Cartesian 
coordinates according to the definition by Setyawan 
and Curtarolo 0. 


A representation in Cartesian coordinates of the 
“standard” conventional cell, excluding hexagonal 
{hP) and rhombohedral [hR) cells, is 

a'3^ = (a',0,0), 

a' 2 ^ = ( 6 ' cos 7 ', b' sin 7 ', 0 ), 

^3 
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where 


= c' cos /3, 

, (cos a! 
= c- 


— cos /3' cos 7 ') 
sin 7 ' 


c 


/ 

2 



— cos^ a' — cos^ /?' + 2 cos a' cos j3' cos 7 ' 
sin 7 ' 


For hexagonal cells, the representation is 



a'^ = ( 0 , 0 ,c'), 


and for rhombohedral cells, 


/I If ^ ^ A \ 

= [a cos—,-a sin—,0 ) , 

a 2 = ( a cos—,a sin—,0 

,T / ,cosa' , f 

^3 = «-^,0,a ./I - 

\ cos -;r V COS 


2 al 


This means that only three cases are necessary in 
a code to output the “standard” conventional ba¬ 
sis vectors in Cartesian coordinates, and there is 
no need to provide individual cases for all 14 Bra- 
vais lattices. Reduction of cases is available for the 
“standard” primitive cell as outlined in the follow¬ 
ing section. 


5.2. ’’Standard” primitive cell 

The transformation matrix P converting “stan¬ 
dard” conventional to “standard” primitive cells is 
defined as 


(api, ap2, apg) — (a',, a^, a3)P. 

The “standard” primitive cell is the same as the 
“standard” conventional system for Bravais lattices 
without centring, and therefore P is the identity 
matrix. 


5.2.1. mS 

For monoclinic base-centred lattices {mS, note 
the definition is different from oS), 


a 


PI 

, T 
P 2 




2 ’ 2 


ap 3 = ( 0 , c' cos a', c' sin a '), 

1 1 /l 1 0 ^ 


P = 



p-^ = 



( 10 ) 


5.2.2. oS 

For orthorhombic base-centred lattices {oS, note 
the definition is different from mS), 

ipi 

T ^ 

^P 2 


2 ’ 2 


' 2 ’ 2 ’^^ 


^P 3 “ ( 0 ; 0 ; C ) , 



(h h o\ 


ll—1 

1—I 

p = 

1 1 0 


110 


Vo 0 l) 


lo 0 l) 


( 11 ) 


5.2.3. ol, tl, and cl 
For body-centred lattices (oJ, tl, cl), 

, T_f a' b' c'\ 

2 ’ 2 ’ 2 y ’ 

, T_fa' _b' c'\ 

V 2 ’ 2 ’ 2 y ’ 

2 ’ 2 ’ 2 y ’ 


p = 

(1 1 i\ 

2 2 2 

1 i 1 

2 2 2 

, P- = 

/o 1 l\ 
10 1 


U 1 iJ 
\2 2 2/ 


1 o; 


5.2.^. oF and cF 
For face-centred lattices {oF, cF), 



( 12 ) 
















p = 


/o 


1 

2 



1 

2 

0 

1 

2 



(1 1 i\ 

1 T 1 . (13) 

VI 1 V 


6. Lattice constant dependence on BZ topol¬ 
ogy 

The topology of the BZ depends on the lat¬ 
tice constants in some Bravais lattices. The con¬ 
ditions that decide the topology, using crystallo¬ 
graphic conventional and “standard” Q conven¬ 
tional lattice constants, are outlined below. Special 
cases for “virtual” crystals are omitted. 

6.1. cP, cF, cl, tP, oP, ol, oS, hP, and mP 
There is only one topology for each Bravais lat¬ 
tice. 


6.2. tl 

There are two topologies for the body-centred 
tetragonal lattice as given in Table [51 

6.3. oF 

There are two topologies for the face-centred or¬ 
thorhombic lattice. The condition for the crys¬ 
tallographic conventional cell depends on the or¬ 
der of basis vector lengths and is rather involved. 
The “standard” conventional cell relabels the crys¬ 
tallographic conventional basis vectors such that 
a' < b' < c', and the condition for each topology 
can be written easily using “standard” conventional 
basis vector lengths as given in Table [HI Then the 
condition for the crystallographic conventional cell 
is obtained using that for the “standard” conven¬ 
tional cell and Table |3l 


6 . 4 . hR 

There are two topologies for the rhombohedral 
lattice as given in Table (TUI 


6.5. mS 

There are three topologies for the base-centred 
monoclinic lattice as given in Table [TTl The crystal¬ 
lographic conventional cell is C-centred with unique 
axis b and /? > 90°, while a “standard” conventional 
cell is C-centred with unique axis a and a' < 90°. 
The quantity that appears in the conditions can 
be expressed using “standard” conventional lattice 
constants as 


cos k'p^ = 


a/2 _ y2 

q /2 _|_ 5/2 


The reasoning behind how the conditions are deter¬ 
mined is provided in [Appendix B| 


6.6. aP 


The reciprocal “standard” conventional cell for a 
triclinic lattice is a Niggli reduced cell in reciprocal 
space where the axes are relabeled such that the in- 
teraxial angles are either all larger (all-obtuse) than 
or all smaller (all-acute) than 90° and is the clos¬ 
est to 90°. Nevertheless there is only one topology 
of the BZ, that is, the truncated octahedron. The 
proof of this is given in [Appendix C[ There is a 
distinct difference between the all-obtuse and all¬ 
acute reciprocal “standard” conventional cells re¬ 
garding the orientation of reciprocal basis vectors 
in the BZ. All reciprocal basis vectors penetrate 
hexagonal faces in an all-obtuse reciprocal “stan¬ 
dard” conventional cells. Our investigation of BZs 
of all-acute triclinic structures using the Materi¬ 
als Project database [l^ shows discrepancies from 
the claims that Setyawan and Curtarolo Q make. 
Setyawan and Curtarolo writes that the BZ is a 
truncated octahedron and the faces of the BZ where 
“standard” reciprocal basis vectors a'j*, a^*, and a'^* 
penetrate are hexagonal, hexagonal, and quadrilat¬ 
eral (parallelogram), respectively. We find that a!* 
penetrates a parallelogram face and a!.^ and ag* pen¬ 
etrates hexagonal faces in the structure of Ag 2 S, 
k'^ = 0.221 A-i, fc' = 0.229 A-i, fc' = 0.105 A-^, 
k'^ = 80.80°, = 79.36°, and = 83.58°. On the 

other hand, a^* penetrates a parallelogram face and 
a'j* and ag* penetrates hexagonal faces in the struc¬ 
ture of HgOsCuSe, < = 0.184 A-^, k[ = 0.207 
A-i, fc' = 0.141 A-i, k'^ = 71.00°, k '0 = 75.21°, 
and /c(, = 77.19°, while ag* penetrates a parallel¬ 
ogram face and a(* and a^* penetrates hexagonal 
faces in the structure of B 307 Rb 2 , k'^^ = 0.095 A“^, 
fc; = 0.102 A-\ k'^ = 0.156 A-i, k'^ = 78.58°, 
k'p = 76.56°, and = 88.15°. In short, always 
one of {a(*, a^*, ag*} penetrates a parallelogram face 
and two penetrates hexagonal faces when the “stan¬ 
dard” cell is all-acute, but which one penetrates a 
parallelogram depends on the lattice constants. In 
light of this complexity, it is recommended to use 
the parallelepiped of reciprocal basis vectors [1, @| 
to describe band diagrams of triclinic structures un¬ 
less there is a very strong reason not to do so. 
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Table 8: BZ topology of a body-centred tetragonal lattice. 


Crystallography 

“Standard” 

Topology 

c > a 

c' > a' 

Truncated octahedron 

c < a 

c' < a' 

Elongated dodecahedron 


Table 9: BZ topology of a face-centred orthorhombic lattice. 


Crystallography 

“Standard” 

Topology 

See the text in Sec. 16.31 

1 ,1,1 

Truncated octahedron 

See the text in Sec. 16.31 

1 ^1,1 
^ ^ ^ -1- 

Elongated dodecahedron 


7. Issues regarding special k-vector point 
definition 

The philosophy behind labeling of reciprocal 
space is very different between the database on 
the Bilbao Crystallographic Server i i i and 
Setyawan and Curtarolo Q. In short, the Bilbao 
Crystallographic Server defines and uses a set of la¬ 
bels for each type of representation domain, which 
is described afterwards, and BZ topology. How¬ 
ever, fractional coordinates of irrational special k- 
vector points are not provided. On the other hand, 
Setyawan and Curtarolo provides a set of special 
k-vector point symbols, including irrational special 
k-vector points, for every BZ topology in each Bra- 
vais lattice. The drawbacks are that their “stan¬ 
dard” conventional basis vectors are different from 
the crystallographic conventional basis vectors in 
many cases and that their definitions are not based 
on the representation domain. 

A number of concepts have to be introduced 
before looking into how labels of reciprocal space 
are defined in the Bilbao Crystallographic Server 
SSi- The reciprocal lattice is represented us¬ 
ing basis vectors in reciprocal space, {a*,a 2 ,a 3 }. 
The reciprocal space group, which is defined as a 
semidirect product of the point group and the trans¬ 
lational group of the reciprocal lattice, is isomor¬ 
phic with the symmorphic space group. This means 
that discussion on a reciprocal space group type can 
be carried out using the corresponding symmorphic 
space group type in direct space. As a result, con¬ 
cepts defined in direct space, for instance Wyckoff 
positions and the asymmetric unit that are defined 
in the ITA Q , can be applied to the reciprocal space 
group. A representation domain, which is a sim¬ 
ply connected parts of the BZ that contains exactly 


one reciprocal space vector (k-vector) of each or¬ 
bit of k, is assigned to each reciprocal space group 
type. The definitions of all 73 representation do¬ 
main types are available on the Bilbao Crystallo¬ 
graphic Server. Orbits of k in reciprocal space can 
be categorized into symmetry-equivalent k vectors 
that correspond to point orbits, or Wyckoff posi¬ 
tions, of the relevant direct space symmorphic space 
group 0. 

Labels of reciprocal space for every representa¬ 
tion domain type are defined in the Bilbao Crys¬ 
tallographic Server i i, i using Wyckoff posi¬ 
tions of the corresponding direct space symmorphic 
space group type. Wyckoff positions can have zero, 
one, two, or three independent coordinate variables; 
therefore, labels are defined for special k vector 
points, lines (line segments, in practice), planes, or 
general positions (GP), respectively. The number 
of labels per Wyckoff position is not necessarily one. 
For instance, coordinates of Wyckoff position la 
(point symmetry mm2) of space group type Pmm2 
(No. 25) can be written as (0,0, z) using coordi¬ 
nates based on the ITA description. Labels dehned 
for this Wyckoff position in the Bilbao Crystallo¬ 
graphic Server are r(G'M) {z = 0), Z {z = 1/2), 
k{LD){Q < z < 1/2), and LE (-1/2 < z < 0) 
iii. The names of labels depend on the BZ 
topology in body-centred tetragonal, face-centred 
orthorhombic, and rhombohedral Bravais lattices 
and the relation between basis vector lengths in 
base-, body-, and face-centred orthorhombic Bra¬ 
vais lattices. The topology of the BZ and relative 
orientation of basis vectors are not explicitly con¬ 
sidered in monoclinic, base-centred monoclinic, and 
triclinic Bravais lattices, and therefore only one set 
of labels are defined per representation domain. 

Band diagrams are drawn along a path in recipro- 
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Table 10: BZ topology of a rhombohedral lattice. 


Crystallography 

“standard” 

Topology 

\/3a < \/2c 

O 

O 

V 

Truncated octahedron 

VSa > \/2c 

a' > 90° 

Rhombic dodecahedron 


Table 11: BZ topology of a monoclinic lattice. 


Crystallography 


b < a sin P 
h > asin/3, 
acos/3 


aP sin^ I 


c 6^ 

b > a sin /3, 

a cos P sin^ 


J£l. 


< 1 


> 1 


“Standard” 


4 ^ > 90 ° 

< 90°, 

b' cos a' b'"^ sin^ a 


c o! ^ 

< 90°, 

h' cos oP 6'^ sin^ o! 


< 1 


> 1 


Topology 

Truncated octahedron 
Elongated dodecahedron 


Truncated octahedron 


cal space. The path is typically along line segments 
connecting special k-vector points on the BZ sur¬ 
face as well as line segments connecting the T point 
and special k-vector points. As a result, symbols 
of special k-vector points on the BZ surface are the 
only relevant labels of interest besides the T point 
in band diagram calculations. In addition to la¬ 
bels, information is also necessary on fractional co¬ 
ordinates of special k-vector points where the band 
diagram path bends. 

The fractional coordinates of centers of faces are 
always rational in a “standard” cell because the k- 
vector of a face center has the form 

|api + |ap 2 + ^a'p *3 where p,( 7 ,r e {-1,0,1}. 

In contrast, some BZ vertices, and therefore cen¬ 
ter of edges, are irrational special k-vector points 
in a number of Bravais lattices (see the definition 
of irrational special k-vector points at the end of 
Section [2]). Table fT^ shows the distribution of Bra¬ 
vais lattices and BZ types among the 58,055 crystal 
structures in the Materials Project database [l^ 
(obtained May 25, 2015). We hnd that crystals 
with irrational BZ vertices account for 59.8% of 
the structures in this database. Still 46.6% remains 
even if we choose to not consider triclinic cells. 

One might argue that forcing usage of the first 
BZ is not always a good idea when drawing a band 
diagram because line segments connecting special 
k-vector points are not necessarily on high symme¬ 
try lines and, therefore, band paths should span 
the asymmetric unit that is uniquely defined for 


each representative domain. However, taking this 
viewpoint does not mean that defining symbols and 
identifying fractional coordinates of irrational spe¬ 
cial k-vector points is unnecessary. For example, 
consider space group type 74 (No. 79) with c/a < 1 
(elongated dodecahedron BZ). The definitions of la¬ 
bels in the Bilbao Crystallographic Server [1, d, @| 
and coordinates according to the ITA description 
is used in this paragraph. Line segment A{LD) 
consists of points with coordinates (0,0, z) where 
0 < z < Mq. The picture of the BZ on the Bil¬ 
bao Crystallographic Server shows that line seg¬ 
ment A{LD) starts at the P [GM) point and ends at 
point Z, and points with coordinates (0, 0, z) with 
Ido ^ z <1/2 form line segment VAi that connects 
points Z and Mq. The quantity Ido is very im¬ 
portant. Line segments A{LD) and VAi together 
form a straight line segment in reciprocal space and 
come in touch at point Z at coordinate (0, 0,Wo)- 
However, the two line segments differ significantly 
in character because line segment A{LD) connects 
the P {GAI) point and point Z at the surface of 
the BZ while line segment VAi is entirely on a BZ 
edge. The quantity Ido and the definitions of points 
Z and Mq as well as line segment VAi are not 
provided in tables on the Bilbao Crystallographic 
Server. Although we can deduce using Setyawan 
and Curtarolo Isl that 


unfortunately, cumbersome conversion of basis vec- 
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tors and coordinates are necessary between the Bil¬ 
bao Crystallographic Server and Setyawan and Cur- 
tarolo. The coordinates of special k-vector points 
in crystallographic conventional cells may be de¬ 
rived from those of “standard” primitive cells us¬ 
ing transformation matrices. Namely, if “standard” 
primitive basis vectors are related to those of crys¬ 
tallographic conventional cells through 

(api)ap2:ap3) = (a'l. ^2, aa)^ = (ai, a2, a3)MP, 

(14) 

the basis vectors of the reciprocal crystallographic 
conventional cell are related to those of the recip¬ 
rocal “standard” primitive cell by 


/aA 


^api\ 

a2 

= MP 

ap2 





Naively mixing symbols in the Bilbao Crystal¬ 
lographic Server M i i and those in Setyawan 
and Curtarolo Q is out of the question. Exist¬ 
ing symbols attributed to the same special k-vector 
point may or may not differ based on the two def¬ 
initions; therefore, the definition must be clearly 
stated when describing results. The symbol for the 
same point differ for example in space group type 
R3m (No. 160) with '/Sa < \f2c (rhombic dodeca¬ 
hedron BZ). Here, point T in the Bilbao Crystallo¬ 
graphic Server corresponds to point Z in Setyawan 
and Curtarolo, and the label P is used to indicate a 
line segment in the former whereas P is a BZ edge 
center in the latter. Existence of multiple change 
of basis matrices for some Bravais lattices gives rise 
to another problematic situation. Consider the or¬ 
thorhombic space group type Pmmn (No. 59). The 
labels of BZ vertices for the same fractional co¬ 
ordinates are the same between the Bilbao Crys¬ 
tallographic Server and Setyawan and Curtarolo, 
therefore the same special k-vector point symbol 
represents the same point in the two definitions as 
long as the crystallographic and “standard” con¬ 
ventional cells are identical. Basis vector lengths 
satisfy a < 6 < c in VOBr (basis vector lengths: 
a = 3.473 A, 6 = 3.854 A, c = 9.050 A) hence crys¬ 
tallographic and “standard” conventional cells are 
the same and the special k-vector point symbols 
are the same in the two definitions. In contrast, 
in CdCu 202 (basis vector lengths: a = 4.002 A, 
b = 9.889 k, c = 3.680 A) basis vectors are re¬ 
taken in the “standard” conventional cell such that 
a' < b' < c', therefore the same special k-vector 


point symbol specifies different points in the BZ. 

In summary, the following two issues need to be 
addressed when performing band diagram calcula¬ 
tions. First, a unique choice of basis vectors must 
be determined for band diagram calculations. The 
change of basis matrix need to be readily available 
if the choice of basis vectors to be used in band 
diagram calculations is different from the crystallo¬ 
graphic convention. Second, a suggested path to¬ 
gether with definitions and fractional coordinates of 
special k-vector points on the path should be prede¬ 
termined for each representation domain type and 
BZ topology. Time-reversal symmetry, in which 
band positions at k becomes the same as at —k, 
significantly reduces the number of representative 
domain types that have to be considered. Time- 
reversal symmetry forces inversion symmetry to ex¬ 
ist in the representative domain type of the BZ re¬ 
gardless of whether the original representative do¬ 
main type has inversion symmetry or not. A list of 
symmorphic space group types with inversion sym¬ 
metry is given in Table [T51 In other words, if there 
is time-reversal symmetry the suggested band di¬ 
agram paths in reciprocal space need to be con¬ 
sidered only for representative domain types corre¬ 
sponding to the 24 space group types in Table [13] 
instead of all 73 representative domain types. 

8. Summary 

An algorithm to obtain “standard” primitive cells 
for efficient and systematic band diagram calcula¬ 
tions is derived. First, a crystallographic conven¬ 
tional cell based on the definition outlined in this 
work is obtained, for example, by using the spglib 
code [l^ versions 1.7.3 and higher. Next, atom po¬ 
sitions are obtained in Cartesian coordinates and 
basis vectors are retaken as necessary to derive the 
“standard” conventional cell. Finally, the atom po¬ 
sitions are converted to fractional coordinates and 
the basis vectors of the “standard” primitive cell 
are determined in Cartesian coordinates. Band di¬ 
agrams often pass through special k-vector points 
with irrational coordinates, and we find that there 
is a need for a crystallographic convention-friendly 
scheme to describe such irrational special k-vector 
points. 
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Appendix A. Derivation of M for mS 

We provide simple proof that one change of basis 
matrix suffices for all base-centred monoclinic cells. 
The reciprocal primitive basis vectors of a base- 
centred monoclinic lattice can be taken in Gartesian 
coordinates as 


ap*i = (4',<?,0), 

ap2 = 

= ( 0 ,-(r-f n( 3 '),s) 


using four positive variables p, q, r, and s and an 
integer n. The restriction that this reciprocal cell 
is primitive uniquely determines p, q, and s. The 
quantity r can be unambiguously defined such that 
r < q, and the cell is “standard”, or ap- = ap*, 
when n = 0. The direct space basis vectors of the 
corresponding primitive cell can be expressed as 


ipi, ap2. 




1 

2 



r-\-nq 

qs 



The basis vector lengths and interaxial angles of the 
corresponding conventional cell with basis vectors 
(af ,ar,a"*) are 


is used to transform this a-axis, C-centred cell 
to a 6 -axis unique, C-centred cell with ba¬ 
sis vectors (a"', a^", ag") through (a",a 2 ,a 3 ) = 
{a'C, 3 . 2 'The basis vectors are minimized 
in a crystallographic conventional cell, therefore 
n = 0. This is exactly the condition for a “stan¬ 
dard” conventional cell. 


Appendix B. BZ topology of mS 


We present the reasoning behind how conditions 
to determine the topology of the BZ of base-centred 
monoclinic cells are determined. This section uses 
the shorthand notation b' = ap^ and discussion is 
based on “standard” primitive cells. 

The reason why the shape of the first BZ of base- 
centred monoclinic cells can be categorized into 
three types is discussed. Gases where inequality 
relations become exact equalities are not consid¬ 
ered. We rely on the fact that the possible shapes 
of Voronoi cells are limited 17| and the following 
lemma: 

Take two reciprocal “standard” primitive basis 
vectors b' and b', where i,j G {1,2,3} and i ^ j. 
The two-dimensional (2D) BZ in the plane where 
b' and b' lies is a hexagon and must: 
a) if b' • b' > 0 , the points 


b' b' b'-b' 

T’Y’ 2 


and 


^ Vir + nq)"^ + 
qs 

c" = 1, 

s 

„ // + m) 

IS a = — , , 

Y(7' + nqY + s^ 

P" = 90°, 

7 " = 90°. 


The set of basis vector lengths is invariant when the 
change of basis matrix 


M = 


/O T 
1 0 
\0 0 


and these negatives exist on centers of different 
sides of the 2DBZ, or 
b) if b' • b' < 0 , the points 

b' b' b' -f b' 

Y’Y’ 2 

and these negatives exist on centers of different 
sides of the 2DBZ. 

As b} • bg = b 2 • bg < 0 by definition, 

b; b '2 b'g b}-fb(, h '2 + h'^ 

2 ’ 2 ’ 2 ’ 2 ’ 2 

and these negatives automatically become centers 
of faces. This means that the BZ must have 10 
or more faces, which forces the BZ to be either an 
elongated dodecahedron (faces are four hexagons 
and eight parallelepipeds), a rhombic dodecahedron 
(faces are 12 parallelepipeds), or a truncated oc¬ 
tahedron (faces are eight hexagons and six paral¬ 
lelepipeds) . 
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1) If b'^ • b 2 < 0, or kp^ > 90°, then and its 

negative will be centers of faces. This is geomet¬ 
rically possible only when the BZ is a truncated 
octahedron and ^ and these negatives are 

at the center of hexagonal faces. and its 

negative are also on a hexagonal face. 

2) If b'^ • b 2 > 0, or kp^ < 90°, there are two cases. 

One is when jg g^gj ^j^g g2 

is an elongated dodecahedron, and the remaining 
case is when is on a face and the BZ is a 

truncated octahedron. 


Appendix C. BZ topology of aP 


We show that the BZ must be a truncated octa¬ 
hedron if there is no pair of basis vectors that are 
perpendicular to each other. According to Horvath 
[l7|, the five topologies that a BZ can take is the 
parallelepiped, hexagonal prism, rhombic dodeca¬ 
hedron, elongated dodecahedron, and truncated oc¬ 
tahedron. The concept of a “zone” is introduced, 
where any edge e of a Brillouin zone determines a 
zone of faces in which each face has two sides equal 
and parallel to the given edge e. The number of the 
opposite pairs of relevant faces corresponding to a 
given zone is two or three, and if it is two then the 
corresponding faces are orthogonal to each other 
This means that a pair of basis vectors of the 
Brillouin zone must be perpendicular to each other 
if the number of pairs of faces in a zone is two. In a 
triclinic cell, basis vectors are generally not perpen¬ 
dicular to each other and therefore all zones should 
consist of three pairs of faces. Edges that are shared 
by hexagons in an elongated dodecahedron consti¬ 
tute a zone of two pairs of hexagons, thus the BZ 
of a triclinic cell cannot be an elongated dodecahe¬ 
dron. The rhombic dodecahedron is a special case 
of the elongated dodecahedron where the length of 
edges shared by hexagons is 0, and at least one rect¬ 
angular cross-section where the centers of edges are 
penetrated by basis vectors always exist. Therefore, 
the rhombic dodecahedron also cannot be a BZ of 
a triclinic lattice. One can easily show that zones 
with two pairs of faces exist in a hexagonal prism 
and in a parallelepiped, hence these also cannot be 
BZs of a triclinic cell. In summary, the BZ of a 
triclinic cell must be a truncated dodecahedron. 
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Table 12: Breakdown of 58,055 crystals in the Materials Project by Bravais lattice and BZ topology. Symbols 
denoted with an asterisk (+) indicate Bravais lattices where some BZ vertices are irrational special k- 
vector points. BZ topologies are TO: truncated octahedron, ED: elongated dodecahedron, and RD: rhombic 
dodecahedron. All-obtuse and all-acute cells are distinguished for triclinic cells. “Other” in base-centred 
monoclinic and triclinic cells indicate that an reciprocal interaxial angle which should be larger or smaller 
than 90° is very close to 90° and therefore the cell cannot be categorized into any case. 


Crystal system 

Centring 

Symbol 

Topology 

Count 

Ratio 

Triclinic 

Primitive 

aP* 

All 

7646 

13.2% 




(All-obtuse) 

3831 

6 .6% 




(All-acute) 

3797 

6.5% 




(Other) 

18 

0 .0% 

Monoclinic 

Primitive 

mP* 


8858 

15.3% 


Base-centred 

mS* 

All 

6919 

11.9% 




b < a sin /3 TO 

4973 

8 .6% 




b > a sin P ED 

1288 

2 .2% 




b > a sin /3 TO 

641 

1 .1% 




Other 

17 

0 .0% 

Orthorhombic 

Primitive 

oP 


7572 

13.0% 


Base-centred 

oS* 


2731 

4.7% 


Face-centred 

oF* 

All 

476 

0 .8% 




ED 

338 

0 .6% 




TO 

138 

0 .2% 


Body-centred 

ol* 


1080 

1.9% 

Tetragonal 

Primitive 

tP 


2696 

4.6% 


Body-centred 

tl* 

All 

3544 

6 .1% 




ED 

826 

1.4% 




TO 

2718 

4.7% 

Hexagonal 

Primitive 

hP 


5496 

9.5% 

Rhombohedral 

Primitive 

hR* 

All 

3435 

5.9% 




TO 

2861 

4.9% 




RD 

574 

1 .0% 

Cubic 

Primitive 

cP 


2096 

3.6% 


Face-centred 

cF 


4775 

8 .2% 


Body-centred 

cl 


731 

1.3% 
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Table 13: List of symmorphic space group types with inversion symmetry. Brackets indicate the space 
group number. 


Crystal system 

Centring 

Symbol 

Space group type 

Triclinic 

Primitive 

aP 

Pi (2) 

Monoclinic 

Primitive 

mP 

P2/m (10) 


Base-centred 

mS 

C2/m (12) 

Orthorhombic 

Primitive 

oP 

Pmmm (47) 


Base-centred 

oS 

Cmmm (65) 


Body-centred 

ol 

Immm (71) 


Face-centred 

oF 

Fmmm (69) 

Tetragonal 

Primitive 

tP 

PA/m (83), PA/mmm (123) 


Body-centred 

tl 

/4/m (87), lA/mmm (139) 

Hexagonal 

Primitive 

hP 

P3 (147), P31m (162), P3ml (164), 
PQ/m (175), P6/mmm (191) 

Rhombohedral 

* 

hR 

P3 (148), P3m (166) 

Cubic 

Primitive 

cP 

Pm3 (200), Pm3m (221) 


Body-centred 

cl 

I m3 (204), Im3m (229) 


Face-centred 

cF 

Pm3 (202), Fm3m (225) 


* Primitive or triple hexagonal. 
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